DVCS and the skewness effect at small x 
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We analyze small- a; DVCS data using flexible GPD models and compare our outcome 
with the full Shuvaev transformation. We point out that the full Shuvaev transform is 
a model that is equivalent to a conformal GPD and a minimalist "dual" parameteriza- 
tion. Some mathematical subtleties of conformal representations are recalled. We also 
comment on a speculation of a factorization breakdown in DVCS. 



Deeply virtual Compton scattering (DVCS), 7* (qi) p(P\) — ► 7(92) jK^), is viewed as the 
cleanest process to access generalized parton distributions (GPDs), which encode a partonic 
description of the nucleon structure. In the kinematics of HI and ZEUS collider experiments 
at HERA, the DVCS cross section is to a large extent dominated by the flavor singlet part 



^ r~ I of the helicity conserved Compton form factor (CFF), denoted as ^i.: 

%(WAQ 2 )*^fe\ s n^t = A*,Q*)\ 2 . (i) 
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Here, W is the c.o.m. energy, A = P 2 — P\ is the momentum transfer, — Q 2 = q\ is the 
incident photon virtuality, and £ = Q 2 / (2W 2 + Q 2 ) « XBj/2 is a Bjorken-like scaling variable. 

The CFF °fl factorizes further into a convolution of the partonic, i.e., hard scattering, 
amplitude C = ( S C, G C) and GPDs H = (^H, G H), (S=singlet quark, G=gluon), 

O : r i 

ON. S H(ti,t,Q 2 )= / dxC(x,tQ 2 /v 2 ,aM)H(x 1 r 1 = Z,t 1 iJ 2 ) , (2) 
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where the skewness parameter 77 = —A • q/ (Pi + P2) ■ q is set equal to £. The factorization 
scale [i separates short- and long-distance dynamics and is often taken as \i = Q. The scale 
' dependence is governed by evolution equations. Note that gluons do not directly enter the 

DVCS amplitude at leading order (LO), but rather drive the evolution of singlet quarks. 
Since the momentum fraction x is integrated out in @, GPDs cannot be directly revealed. 

It is our objective to find flexible GPD models, which satisfy the known theoretical 
constraints, and which can be pinned down by fits to HI and ZEUS DVCS data [2] at LO 
and beyond [5]. Thereby, we find it convenient to work with conformal GPD moments. For 
integral conformal spin j + 2 they are defined by convolution with Gcgenbaucr polynomials 
Cj(x), e.g., for quarks: 
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where q is the flavor index. In the forward limit A — > these moments simply reduce for 
j = 0, 1,2, • • ■ to familiar qj(^ 2 ) Mcllin moments of parton distribution functions (PDFs). 
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Conformal symmetry guarantees that they evolve autonomously under evolution at LO (ex- 
cept for the quark-gluon mixing) . Now we require an appropriate behavior of the conformal 
moments ([3]) for j — ► oo with |arg(j)| < ir/2 [3]. Carlson's theorem then states that their 
analytic continuation with respect to j is unique and the GPD moments can be used to 
calculate the corresponding CFF, cf. 
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C J (Q 2 /M 2 ,«s(M))^ J (e,t,M 2 ) (4) 



within a Mcllin-Barnes integral (3j HJ \5\ . Here the singularities of the integrand, except for 
those of tan(7rj/2) at j = 1, 3, • • • , lie on the l.h.s. of the integration path. 

To get flexible GPDs and CFF H and to make (loose) contact with Regge terminology 
(UH], we expand their moments in terms of t-channel SO (3) partial waves (PWs) dg [6], 
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(5) 



which are labelled by the i-channel angular momentum J. These PWs can be expressed 
in terms of the familiar Legendre polynomials, where the cosine of the scattering angle 9 
may be approximated by —1/rj and which are normalized to one for i] — > 0. We include an 
effective "pomcron" pole at a(t) = ao + a't in the PW amplitudes and parameterize the 
residual t-dependence f3j(t) by a dipole or exponential ansatz. The strengths hj of leading 
J = j + 1 PWs and the intercept a = a(t = 0) are obtained from fits to some DIS data, 
and the remaining parameter space is further reduced by either truncation of the sum §E§ or 
its model dependent resummation, so-called S-PW model, and is constrained by DVCS fits, 
which turn out to have good quality % 2 /d.o.f. w 1, see Fig.[T] (We take jtio = 2 GeV as input 
scale.) One important partonic quantity that can be extracted from such fits is the skewness 




Figure 1: LO fit to the DVCS cross section from HI and ZEUS versus W (left) and Q 2 (right), 
ratio of GPDs, which is given by the value of the skewness function at = r\Jx = f : 

r(x, tf, /i ) = — — — r—oT-. (6) 



Before we present results of our fits, some issues should be clarified. Let us recall that 
there are at least four different conformal GPD representations that are general and are 
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used in phenomenology: Shuvaev-NoritzscfQ transformation [7J [5] , "dual" parameterization 
[5], and two versions of Mellin-Barnes representations [31 [TU]. All of them start from an 
integral conformal SL(2 |R) PW expansion of the (crossed) GPD in terms of Gegenbauer 
polynomials (now viewed as generalized functions in the mathematical sense), 



H*(x, v ,t) 
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however, in the end they provide different representations for GPDs. Since these specific 
conformal representations involve some mathematical subtleties, efforts to understand them 
have been undertaken in Refs. [HI [21 [TO1 HU E]. We consider the one-to-one correspon- 
dence among the various representations established already by construction, as sketched 
by the thick arrows in Fig. [2l although the in- 
verse transformations are only partly known. 
Furthermore, if some mathematical subtleties 
are ignored, general conformal representations 
degenerate to parameter-free GPD models for the 
small-x region. 

In the small- x region the one-to-one corre- 
spondence among the different parameter-free 
models can be analytically shown by noting that 
they all have the conformal value of skewness 
function, e.g., for t — and fixed /j, 2 
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Figure 2: The one-to-one relation of various 
conformal GPD representations is indicated 
by thick arrows. The restriction of minimal- 
ist GPD models is sketched by thin arrows. 



Assuming, like in most standard PDF paramc- 
tcrizations, that effective Rcggc behavior is in- 
cluded in the intercept a = a(t = 0), one 
easily finds the skewness function ([H]) from the 
Mellin-Barnes representation [3j [10] of a confor- 
mal GPD model [4] by shifting the original integration path to the l.h.s. in the complex 
j-plane, cf. Eq. ([4]). Thereby, the hypergeometric function in ([8]) is nothing but the Clebsch- 
Gordan coefficient of the conformal PW expansion, taken at j = a — 1. As long as the 
conformal moments behave smoothly in the vicinity of rj = 0, a GPD model degenerates to 
a conformal one. For the full Shuvaev transformation model the skewness function ([5]) can 
be found in Ref. |13j as the standard integral representation of hypergeometric functions, 



2 2a+1 r(3/2 

r(i/2)r(n 



Jo 



(9) 



For an effective "pomeron" (a ~ 1) and "Reggeon" (a ~ 1/2) pole such analytic approx- 
imation works quite well for x — £ < 10~ 2 . Further singularities (0 < a < 1) might be taken 



a This integral transformation was proposed by A. G. Shuvaev in [?]. The term (full) Shuvaev transfor- 
mation was used in Refs. 1121 1131 for an integral transformation that results in a restricted GPD transform. 
This restriction was removed by J. Noritzsch in Ref. [8]. He was also the first who utilized this integral 
transformation to set up flexible GPD models, describing small-x DVCS data at LO. Thus, — not only to 
avoid confusion — we name the general integral transformation Shuvacv-Noritzsch. 
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into account, too. If the scale grows the effective "pomeron" intercept a increases, according 
to the well-known double-log approximation. If "Reggeon" contributions still play a role, 
the analytic approximation of a GPD will of course be a linear combination of conformal 
skewness functions flEJ. 

As the reader realizes, for the same small- x GPD model two different names are used, 
or even three if one includes the minimalist "dual" parameterization. In the remaining part 
of this presentation we simply adopt the terminology from the "dual" parameterization [9] 
and use, instead of the term "parameter-free prediction", the term "minimalist GPD" for 
both our conformal and the full Shuvaev transformation model. Let us now scrutinize some 
statements presented in Ref. [13], and used there to support the problematic idea that the 
minimalist GPD should be considered as more than just a restricted GPD model [14] , 

i. A naive Taylor expansion of integral GPD moments around ry = would imply that the 
Shuvacv-Noritzsch GPD transform reduces to a PDF at small ij. However, it is well known 
that the operations of the small-77 expansion and of an analytic continuation of polynomials 
with respect to their order do not commute. The Sommcrfcld- Watson transformation of a 
series, representing a GPD in terms of integral moments, is only possible if the assumptions 
of the Carlson theorem are satisfied. The theorem then states that the continuation of 
conformal moments is unique. If the integral GPD moments are first reduced to PDF ones, 
the Sommerfeld- Watson transformation is not applicable [3]. That the Shuvaev-Noritzsch 
GPD transform might not reduce to a PDF is known since Refs. [El [8]. The analogue of this 
in the "dual" parameterization framework is the increase of singular behavior of non-leading 
forward- like functions [11] . 

ii. It was said in JT3] that the absence of 
poles in the right half of complex j plane is 
a physical requirement that is sufficient to 
guarantee that the full Shuvaev tranforma- 
tion (minimalist) model is the unique GPD 
model at small 77 to 0(r]) accuracy. We first 
remark that such a requirement would be 
physical only for the angular momentum J. 
Moreover, counterexamples are known, such 
as resummed S-PW model used here and in 
[5] for DVCS fits. A small-?? expansion re- 
veals then that the poles in the right half of 
complex j plane are an artifact of the expan- 
sion itself and, thus, should be considered as 
spurious [S]. 

Hi. It was stated in |13j that up to order 
0(77) the full Shuvaev transformation is compatible with the NLO evolution equation in the 
MS scheme. However, this is not true, as spelled out, e.g., in [4], and demonstrated in Fig. [3] 
where we show the NLO evolution of a GPD on r\ = x both in the MS scheme and within 
the procedure [13] . which actually corresponds to an evolution operator in CS scheme, as 
dot-dashed and dashed curves, respectively. One realizes that the discrepancy increases with 
evolution and with decreasing x. Hence, it cannot be an 0(rj) effect. 

The phenomenological status of a minimalist GPD model was to a large extent discussed 
in Ref. [5]. It is shown in Fig. 0Ji that a flexible model (solid), pinned down by a good 
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Figure 3: GPD on the cross-over line versus Q 2 
in the MS scheme (dot-dashed) and the CS one 
(dashed). The discrepancy is the error if a full 
Shuvaev transform is used in the MS scheme. 
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Figure 4: Sea quark GPDs on the cross-over line obtained from good fits in CS (solid) and MS 
(dot-dot-dashed) schemes are compared with minimalist GPD models used in our bad xVd.o.f. « 3 
LO fit (dotted) and from Ref. [13] (dashed, dot-dashed) at LO (left) and NLO (right). 



X 2 /d.o.f. ~ 1 fit. is almost two times smaller than a minimalist GPD model, shown by 
the other curves, and so the latter model is ruled out at LO. This was also confirmed 
within the "dual" parameterization and standard PDFs [TB]. (Small shape discrepancy of 
our minimalist model (dotted) and those of Ref. [13] (dot-dashed, dashed) comes from our 
neglecting of "Rcggcon" contributions.) At NLO one might conclude from Fig. |4Jd that 
within experimental and theoretical uncertainties the minimalist models [13] (dot-dashed, 
dashed) are compatible with ours (solid, dot-dot-dashed), obtained from data. Still, a 
more detailed view shows that skewness effect, ^-dependence, and scheme convention are 
interrelated. For instance, we found that in MS scheme the E-PW model with dipole t- 
dependence (dot-dot-dashed) is somewhat away from a minimalist one (which provided only 
X 2 /d.o.f. ~ 1.5). On the other hand in the CS scheme at NLO (solid) and NNLO it turned 
out to be close to a minimalist one. We add that the error from the wrong evolution, 
see Fig. [31 is comparable to the PDF uncertainties and that for the kinematical region of 
interest even the forward evolution operator is not stable in perturbation theory. Note that 
the reparameterization effects for our quark GPD models, resulting from switching on NLO 
corrections, are much larger than for PDFs. We view this as a sign that in our modelling 
we have not yet reached control over evolution. 

In conclusion, the full Shuvaev transform is a minimalist GPD model valid at small- a;, 
which is also realized in other versions of conformal representations. So far we see no general 
theoretical arguments that support preference of this model over others. NLO evolution is 
performed in Ref. [13| in such a way that it is inconsistent with the MS scheme. A minimalist 
model might describe DVCS data at NLO, however, not necessarily when precision level is 
reached. Finally, we recall that the skewness effect for gluons is in all popular GPD models 
approximately the same, namely, zero (i.e. r G (i? = 1) m 1; what is often quoted as a 
large skewness effect is 2 Q ~ 1 ) [5]. Hence, the skewness model uncertainty, compared to the 
phcnomcnological PDF uncertainty, might presumably be considered as small at NLO. Fully 
flexible GPD models, allowing the control over evolution, should be invented and confronted 
with DVCS data in further studies. We hope that more DVCS HERA II run data will 
become available and so the statistical errors will be reduced. This might then allow to 
address the observables and partonic quantities, as discussed in Ref. [5]. 
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Comment added 

At the DIS 2009 conference Felipe J. Llanes-Estrada presented, as collaborator of S.J. Brod- 
sky, J. Tim Londergan, and Adam P. Szczepaniak, a talk |17j with the title "Regge behavior 
in DVCS: non-factorizability and J = fixed pole" in which it was stated that Regge be- 
havior is responsible for the breakdown of collinear factorization. D.M. likes to take the 
opportunity here to repeat his comment, given after the DIS 2009 summary talk "Structure 
Functions and Low-x ", that this speculation is based on a "bad" GPD model. Further clar- 
ifications, e.g., that the illness of this "bad" GPD model arises from the ultraviolet/infrared 
region rather than from the Regge behavior, how to formulate a "good" GPD model in 
the reggeized parton model framework, and how to fix a "bad" GPD model within analytic 
regularization, are presented in Ref. |18| . Interestingly, the authors find it not legitimate to 
use analytic regularization to fix their model, however, suggest the same regularization to 
extract a so-called J = pole within the GPD formalism that, according to them, cannot 
be applied to DVCS (i.e., small —t). In contrast to the collinear factorization approach 
[5], it could not be demonstrated so far that the Regge exchange amplitude model, intro- 
duced two years ago, describes DVCS data. In our opinion the reason is simple, namely, 
the Regge exchange amplitude model is already ruled out by HI and ZEUS data [18] and 
maybe also by the fact that the scaling hypothesis works rather well in a global fit to fixed 
target DVCS measurements [5]. Thereby, a 'partonic dispersion relation' approach is uti- 
lized. Since some clarification is still needed, see conclusion in Ref. [19], we like to spell 
out again, see Refs. [2"0l H] [TS] and references therein, that our uses of 'partonic dispersion 
relations' is in one-to-one correspondence to the collinear factorization approach, taken in 
the LO approximation and supplemented with the scaling hypothesis. 
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